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Abstract 

Primordial magnetic fields may account for all or part of the fields observed 
in galaxies. We consider the evolution of the magnetic fields created by pseu- 
doscalar effects in the early universe. Such processes can create force-free 
fields of maximal helicity; we show that for such a field magnetic energy in- 
verse cascades to larger scales than it would have solely by flux freezing and 
cosmic expansion. For fields generated at the electroweak phase transition, 
we find that the predicted wavelength today can in principle be as large as 
~ 10 kpc, and the field strength can be as large as ~ 1CP 10 G. 



The origin of galactic and intergalactic magnetic fields is an unsolved problem |T|. The 
standard a-Q dynamo theory of galactic magnetic fields || has been criticized for not ad- 
equately taking into account the back reaction of the growing magnetic field, and in any 
case, the theory requires a seed field of unknown origin. There is now an extensive literature 
examining the possibility that magnetic fields were created by exotic processes in the early 
universe, although most such processes, if they work at all, produce fields on scales too small 
to be of interest to astronomy ||. 

In magnetohydrodynamics, energy can be transferred from small to large scales by a pro- 
cess known as the inverse cascade. As shown in pioneering work by Pouquet and collabora- 
tors ||], a critical ingredient of the inverse cascade mechanism is the presence of substantial 
magnetic helicity, also known as the Chern-Simons term (although non-helical cascades have 
also been investigated ||). The idea that magnetic helicity may drive an inverse cascade 
from microphysical magnetic fields to large-scale cosmological fields has been advocated by 
Cornwall ||, and investigated subsequently by Son 0, who proposed scaling properties we 
will verify below. The present work has two goals: 1.) To work towards an analytic under- 
standing of the cascade process as an initial-value problem appropriate for cosmology, by 
studying similarity solutions of the MHD equations in the presence of helicity; and 2.) To 
apply this understanding to proposed mechanisms which create helical primordial fields via 
pseudoscalar processes in the early universe. A preliminary account has been given in |J. 

The early-universe processes we consider can be thought of as arising from the evolution of 
an electrically neutral pseudoscalar field 0, coupled to electromagnetism through a Lagrange 
density 

£=-^0F^ = 0E-B, (1) 

where F^ u = ^e fJ,upa F pa . is the dual electromagnetic field strength tensor. (For simplicity 
we elide the distinction between electromagnetic and hypercharge fields; they are related by 
factors of order unity.) The pseudoscalar in question may represent an axion ||, a more 
general pseudo-Goldstone boson [T^j, or may model the effect of a chemical potential for 
right-handed electron number (in which case the chemical potential fi is given by the time 
derivative of 0). This last scenario has been proposed by Joyce and Shaposhnikov ||TT 



who showed that it could lead to fields at the electroweak phase transition with magnitude 
£>ew ~ 10 22 G and coherent over length scales Aew ~ 10 _8 ^ew ~ 2 x 1CT 9 cm. 

The field equations for electromagnetism in the presence of the interaction (ID) are (in 
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units with c = 1) 

d t E = V x B - J - 0B - V0 x E 

V -E = p E + V0-B 
d t B = -V x E 

V • B = , (2) 

along with Ohm's law, J = cr(E + v x B), where v is the velocity of the fluid. In the 
cosmological context of interest here, is effectively homogeneous (V0 = 0) and the charge 
density pe vanishes. The magnetic field then satisfies 

(V 2 - df)B = a[d t B — V x (v x B)] — V x B . (3) 

For the case that a = it has been shown that solutions of (0) can be unstable to expo- 



nential growth in the magnetic field |TJ]; Garretson, Field and Carroll [HJ considered field 
production by such a mechanism during inflation, concluding that the resulting amplitudes 
were too small to be of astrophysical interest. If one ignores the displacement current <9 t E 
and sets = 0, (0) reduces to the induction equation of dissipative MHD. In the situation 
of interest here the conductivity is non- negligible; the electric field changes slowly on the 
timescales of interest, so we neglect <9 4 E; and the bulk velocity of the fluid is small, so we 
neglect v as well. Hence the appropriate form of (|3[) is 

(d t — t?V 2 )B = r]<j)V x B , (4) 

where rj = l/a is the resistivity. 

It is useful at this point to go to Fourier space and decompose B(k) into modes of definite 
helicity (or equivalently, circular polarization), B(k) = B + u + + B_u_\ here u± — U\ ± iu 2 , 
where iix, «2, and m 3 = k/A; form a right-handed orthonormal basis. Then solutions to 
are of the form 



B±(k,t) = 5 ± (k,0)exp 



-rj J k(k =F 4>) dt 



(5) 



For cf) positive, the B + modes will grow exponentially if k is less than <f), with maximum 
growth rate for k = 0/2, while the 5_ modes decay away. When = the field undergoes 
Ohmic decay, although Joyce and Shaposhnikov show that this effect is unimportant for the 
relevant wavelengths in their scenario [JTTJ . 
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The expectation value of the magnetic energy density of an isotropic plasma in a volume 
V can be expressed in terms of an energy spectrum ej^f as 



E M 



^V- 1 J v B 2 d 3 x = J™ ef dk, (6) 
where 

ef = 2vrfc 2 (B(k)-B*(k)). (7) 

Similarly, the expectation value of the magnetic helicity (or Chern-Simons) density can be 
written 

r roo 
rM _ T r-1 / A T3 ^3 / iM 



H M = V~ / A • B d x = / hfdk , (8) 
Jv J o 

where 

^ = 4vrA; 2 (A(k).B*(k)) . (9) 

Note that while e^ 1 is non-negative, can be of either sign. The helicity and energy spectra 
satisfy a realizability condition, 

K\ < 2k- 1 # . (10) 

We say the field is "maximally helical" if, for every k, is of the same sign and saturates 
this inequality. (See |H| for further discussion.) 



In Coulomb gauge (k ■ A = 0), the modes of the vector potential and the magnetic field 
are related by B± = ±kA±. The energy and helicity spectra are then 



ef = 4nk 2 (\B + \ 2 + \B_ 
hf = 8nk(\B + \ 2 -\B_\ 2 ) . (11) 

Since the B + modes are amplified and the _B_ modes suppressed by the evolution of 0, only 
B + will contribute; such a field satisfies hJjf = 2k~ 1 e^f, and is therefore maximally helical. 

If the spectrum of a maximally helical magnetic field is strongly peaked around some 
wavenumber k p , the configuration will be substantially force-free: J x B = (V x B) x B = 0. 
This can be seen by considering 

V x B(x) = J A;5 + (k)M+e ik x d 3 k ~ fc p B(x) . (12) 
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Force-freedom has been verified in numerical simulations []14"|. This condition plays an 



important role in the evolution of the fields. The velocities associated with force-free fields 
are small, protecting them from the Silk damping (photon viscosity) discussed in [[15] (which 
implies in addition that they will not produce detectable distortions in the spectrum of the 
cosmic microwave background 



We turn now to the principles behind the inverse cascade. Pouquet et al. M studied 
MHD turbulence in the eddy-damped quasi-normal Markovian (EDQNM) approximation, 
in which eddy damping is used to close the nonlinear equations in Fourier space. These 
equations preserve the ideal invariants of total energy and magnetic helicity. Pouquet et 
al. studied the evolution of the spectra e^f, h^, and e\ (the kinetic energy) under various 
conditions. The presence of magnetic helicity causes a shift of magnetic energy from large 
k to small, as the system attempts to minimize magnetic energy while conserving magnetic 
helicity [equation (0)]. In their study of greatest interest here, illustrated in Figure 4 of [§J, 
nonhelical kinetic energy, together with maximally helical magnetic energy, are injected at 
k = 1 at a constant rate, and hJjf is computed for 0.016 < k < 16 at times < t < 200. The 
result is a well-defined wave of magnetic energy and helicity that propagates from k = 1 
downward, reaching k = 0.03 at t = 200. Similar results have been obtained elsewhere 



Here we attempt to quantify this phenomenon analytically so as to generalize it to the 
cosmological case, in which magnetic energy and helicity are injected in a short time. To do 
so, we use equation (3.4) of @j: 

(^Nonlocal = k^T k (hj - k'hf) + f k ef + (a? - af)ef - 2^k 2 hf . (13) 

This equation is obtained by expanding the full EDQNM equations at a given k to isolate 
those nonlinear interactions which are nonlocal in that they involve either k' < ak or k' > 
k/a, where a is a small parameter. The first term in (|I~3"D, with the transport coefficient in 
k - space T k , represents the so-called Alfven effect, in which a relatively strong magnetic 
field at k! << k provides a background on which excitation at k propagates as an Alfven 
wave, bringing into equipartition with e£, and into equipartition with the kinetic 

helicity . The second term is a correction to the Alfven effect. The third term represents 
the a-effect of mean-field dynamo theory, according to which any inequality between the 
kinetic helicity and k 2 h^, the current helicity, at k! » k, tends to destabilize hJjf and e^ 1 
at k. We describe this effect, which is believed to be responsible for stellar magnetic fields 
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T9| , further below. The last term, turbulent viscosity, can be shown to be negligible in our 



case because the velocities are small. 

The third term is of greatest significance here. Particularly important is the appearance 
of the current helicity (b • V x b), whose spectrum is k 2 h]f. The current helicity was 
neglected in early discussions of the a-effect, and was first uncovered by Pouquet et al [|J. 
Its importance has been confirmed in [50] by a very different approach, and it plays a crucial 
role in what follows. 

To proceed, we identify the dominant term in ( |13|) by reference to Figures 4 and 6 of [|J . 
First, note that at the peak k p of the helicity wave at t = 200, we find e[ ~ 0.05e^, as would 
occur if the magnetic field is exerting little force on the fluid. This is consistent with the fact 
that the injected field is maximally helical, hence substantially force-free according to (|12|). 
This view is supported by the fact that in the direct simulation of |14j], the configurations 
formed were still 96% maximally helical (hence substantially force-free) at t = 45, with 
e^/e^ ~ 0.02 at the peak. It is also consistent with the fact that the field at wavenumbers 
less than k p is very weak, and hence unable to bring about equipartition of with e^ 1 
via the Alfven effect .Q We conclude that in our case the Alfven effect is weak, and so we 
neglect the first term in (|T3|), and with it, the next term, a correction to it. An important 
implication is that in the helicity wave, the condition of maximal helicity applies: 

ef = \khf . (14) 

We are left with the a-term. In light of fll4"f) , ( |i~3|) can be written 

(d t hf) NM =~(a y k -af)khf , (15) 

demonstrating explicitly the potential of a\ — a^ 1 to destabilize h^. According to equation 
(3.6) of 1, 

A roo 

«*-«? = -o L BknW - fh?)dq , ( 16 ) 

O Jk/a 



1 The presence of only a very small velocity in Pouquet's solutions raises the question as to how 
the field can evolve in time, absent dissipation. The authors are currently investigating this with 
an approximate model of the Pouquet equation. 
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where the inverse decay constant 9k qq is defined below [equation (0)]. From Figure 6 of 
Q, the lack of kinetic energy in the helicity wave means that \a^\ << \c^f\, and can be 
neglected. Thus 

(d t hf) Nonlocal = —afkhf (17) 

O poo 

= -khf / 6 kqq qX dq . (18) 

O J k/a 

Hence, if is of one sign throughout the wave (as in our case, where we assume that helicity 
of only one sign is injected), an instability is possible. Pouquet et al. (p. 342) describe how 
this instability works. Suppose the wave has reached k p = k\. Because a.k{kv) is due to 
helicity at k > &2, and the helicity is stored at k\, wavenumbers ki < k\ become unstable, 
so that /ijjf and e]f grow there. The growing h]f at k 2 destabilizes k% < k 2 , and so forth. 

This can only be a crude explanation, as it relies on the application of (|I3"D, (|T5[) and (|TBD, 
which are based solely on nonlocal interactions. The reader should bear in mind, however, 
that the numerical results in [|J, represented in Figures 4, 5, and 6 of that paper, are based 
on the full EDQNM equations, not just the nonlocal version in (|i~3|). 

The inverse decay constant 6k qq is given by 

Okqq = GUfc + 2/i,)" 1 , (19) 

where the eddy-damping rate fik is given by 

= iy + v)k 2 + 0.26 ^ p 2 e p dp S 

+ \H k i f* e p d p) ■ ( 2 °) 



1/2 



The first term represents damping by viscosity v and resistivity rj. In astrophysical ap- 
plications v and 7] are typically of order 10~ 6 times the succeeding terms, so this effect 
is important only at very large wavenumbers; we will ignore it henceforth. The second 
term parameterizes damping due to self-distortion, and the third that due to the nonlinear 
interaction of Alfven waves. 

In light of (|I~4"D, we can use ( p3|) to express ( pCD as 

Mfe = (H M e p f 2 F(0 , (21) 

where 
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1/2 

1/2 



F(0 = 0.26 (jf* CMC) d(j 

+0.82£ (j\s(()d(Y . (22) 



We now look for similarity solutions of the form 

hf(t) = g(t)s(0 , (23) 
where the shape function s(£) depends on the wavenumber scaled to the peak: 

* = m- (24) 

We normalize the shape function by J s(£) d£ = 1, which allows us to express the time 
dependence of (f23|) as 

= W ' (28) 

Given this ansatz and H M (t), a solution will be fully specified by the shape s(£) and the 
peak wavenumber k p (t). 

Given (22~), the evolution equation ( |T8| ) can be written 



tfM I ( l\ n x l/2 

- I + = (26) 



# M k p \ ^ s ) V2 



where s' = ds/d^ and 

G «> = HlF(|f§(C) dC - (27) 

Equation fl26|) is valid for ~ k p , in particular for k = k p (^ = 1), where s' = by 
definition. Thus, fc p (t) solves the differential equation 

A / IT M \ / 1 \ 1/2 

^ (4- G(l), (28) 



V /up J v 2 

where G(l) is a dimensionless constant of order unity. Pouquet et al. |4j] considered the case 
H M = const, so H M = H M t, in which case the solution to (B3) is 

k P (t) = k p (U)U/t , (29) 

provided that 



S 



G{1) 



-1/2 



(30) 



The numerical solution of the exact EDQNM equations found in Q fits (29)-(30) very well. 

Reassured that we can find similarity solutions consistent with the numerical results, 
we turn to the H M = case, relevant to cosmology once the pseudoscalar has stopped 
evolving. Then (^) implies that similarity solutions will once again exist, with 

k p (t) = k p (t t ) l^jj , (31) 

provided that 

"9 



G(l) 



-1/2 



(32) 



Note that fl3~T|) verifies a result of Son |7|], derived from a different line of argument. 

In addition to the peak wavelength, we also want to know the rms magnetic field, B rms = 
(B 2 ) 1 / 2 . For our maximally helical fields we have 

(B 2 ) = 2 / ef dk = / khfdk 
Jo Jo 

poo 

= H M k p (t) / is{i) . (33) 
J o 



The integral in the last expression is independent of time. Thus, when H M = constant, ( |31~|) 
implies 

B ims {t) = B rms (U) (J) . (34) 

The preceding discussion has assumed a flat spacetime background; it is straightforward 
to adapt these results to an expanding Robert son- Walker spacetime with metric ds 2 = —dt 2 + 
R 2 (t)dx. 2 = R 2 (t*)[—dt* 2 + dx 2 }, where R is the scale factor and t* is the conformal time. In 
the radiation-dominated era, the complete set of MHD equations is conformally invariant 
(see for example ]2l|]), and our results will apply directly to the conformally transformed 
quantities 

-^rms = R 2 B rras , k* = Rk p . (35) 

Thus, the inverse cascade will be characterized by k*(t*) oc (t*)~ 2 / 3 and B* ms (t*) oc (t*) -1 / 3 , 
where the conformal time in the radiation-dominated era is given by t* = 2(t^Q/ REQjt 1 ^ 2 , 
in which the subscript EQ refers to the epoch of matter-radiation equality. The physical 
quantities at this epoch are therefore related to their initial values by 
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-Brms(^EQ) — -^EQ-^rms (^EQ 
-2 fjl) 



\ r EQ/ 



1/3 
1/6 



and 

M*eq) = -Req^p(^eq) 

/ t * \2/3 



B ims (ti) (36) 



fcp(ti) • (37) 



-Req V^eq/ 

It is possible that additional inverse cascade occurs in the matter-dominated era following 
iEQ- Pending further study, we will assume here that the effect of this is negligible, and so 
the following estimates of the characteristic scale of the field today are lower limits. Thus, 
for t > £eq the field is frozen in, the wavenumber scales as -R -1 , and the magnetic field as 
R~ 2 . Their values today are thus 

1/6 

B rT as(ti) (38) 



and 



kM = ^(±-) 1/3 kM- (39) 

In summary, a maximally helical primordial magnetic field created at time t{ (during ra- 
diation domination), with initial amplitude -B rms (tj) and initial coherence length A(tj) = 
2ir/k p (ti), will undergo an inverse cascade, increasing its length scale by a factor (tEq/U) 1 ^ 3 
over and above stretching due to the expansion of the universe, while its amplitude is diluted 
by an additional factor (ti/tEQ) 1//6 - Note that our estimates of the present coherence length 
substantially exceed those of Son |7[], who terminated the inverse cascade at e + /e~ annihi- 
lation because of increasing viscosity. We have argued that our modes are nearly force-free, 
with negligible velocities; viscosity has no effect. 

For purposes of illustration, let us consider the fate of a magnetic field created at the elec- 
troweak phase transition (T EW = 200 GeV), so tj = £ew = 6 x 10~ 12 sec. We express the ini- 
tial coherence length of the field in terms of the Hubble radius, A(£ew) — /a-^ew = 0-4/a cm > 
and the initial amplitude in terms of the total energy density, -B rm s(^Ew) = /bV^Pew 
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2 x 10 25 /b Gauss, where fx and fs are dimensionless factors less than unity. (We have 
switched here from Lorentz-Heaviside units to CGS in order to express the magnetic field 
in Gauss.) The field today will be coherent over scales 



A(t ) = 5 x 10 22 / A cm = 20/ A kpc , (40) 



with amplitude 



B rms (to) = 4 x W- w f B Gauss . (41) 



In the electroweak case, then, the characteristic length scale of the field has been amplified by 
a factor of (^eqAew) 1 ^ 3 = 6 x 10 7 more than would be expected for a frozen-in configuration. 
If the initial field is coherent over the Hubble radius at the electroweak scale (f\ ~ 1) and 
comparable in energy to the total energy (/g ~ 1), this results in a length scale of ~ 20 kpc 
and an amplitude of ~ 10 -9 Gauss. According to Dolag et al. ||22|| , the primordial field 



required to explain Faraday rotation measures of the Coma cluster is ~ 10~ 9 Gauss. The 
observations are consistent with scales of ~ 60 kpc. Thus, if f\ ~ fs ~ 1, primordial helicity 
could explain the fields in the Coma cluster. 

However, Joyce and Shaposhnikov [p] have estimated that their scenario for helical field 



generation at the electroweak scale results in fields with fx ~ 10~ 8 and fg ~ 2 x 10~ 3 . While 
intriguing (for example as a candidate seed field for a galactic dynamo), these parameters 
fall short of providing a sufficient explanation for the fields seen in galaxies today. It is 
therefore worth considering variations on this mechanism, perhaps with different dynamics 
for the pseudoscalar 0, or field creation at a later epoch such as the QCD scale. Mechanisms 
which create large-amplitude fields without appreciable helicity will undergo significantly 
less (if any) inverse cascade, and thus have a difficult time leading to fields of astrophysical 
significance in the present universe. 



ACKNOWLEDGMENTS 

This work was supported in part by the National Science Foundation under grant 
PHY/94-07195 and NASA under grant NAGW-931, and GF acknowledges discussions at 
the Institute for Theoretical Physics Program on Astrophysical Turbulence, April- June, 
2000. 



11 



[1 



[2 

[3] 
[4] 
[5 



[6] 
[7] 



[9 
[10 

[11 
[12 

[13 
[14 
[15 



REFERENCES 

For reviews see E.G. Zweibel and C. Heiles, Nature 385, 131 (1997); A.V. Olinto, in 
Proceedings of the 3rd RESCEU Symposium, University of Tokyo, Nov. 10-13, 1997, 
Ed. K. Sato, T. Yanagida, and T. Shiromizu; |astro-ph/980705T. 



A. A. Ruzmaikin, A.M. Shukurov, and D.D. Sokoloff, Magnetic Fields of Galaxies, 
Kluwer, Dortrecht (1988). 



O. Tornqvist, |astro-ph/0002307 



A. Pouquet, U. Frisch, and J. Leorat, J. Fluid Mech. 77, 321 (1976). 

A. Brandenburg, K. Enqvist and P. Olesen, Phys. Rev. D 54, 1291 (1996), 
|hep-ph/9608422| ; P. Olesen, |astro-ph/9610154| ; K. Enqvist, |astro-ph/970730q ; T. 



Shiromizu, |astro-ph/9810339 



J.M. Cornwall, Phys. Rev. D 56, 6146 (1997); |hep-th/9704022 



D.T. Son, Phys. Rev. D, 59, 063008 (1999); |hep-ph/9803412 



S.M. Carroll and G.B. Field, in 33rd Rencontres de Moriond: Fundamental Parameters 
in Cosmology, 17-24 January 1998, Les Arcs, France; |astro-ph/9807159. 



M.S. Turner and L.M. Widrow, Phys. Rev. D 37, 2743 (1988). 

W.D. Garretson, G.B. Field, and S.M. Carroll, Phys. Rev. D 46, 5346 (1992); 
[hep-ph/9209238 . 



M. Joyce and M. Shaposhnikov, Phys. Rev. Lett. 79, 1193 (1997); |asFro-ph/9703005 . 



S.M. Carroll, G.B. Field, and R. Jackiw, Phys. Rev. DAI, 1231 (1990); A. A. Andrianov, 
R. Soldati, and L. Sorbo, Phys. Rev. D 59, 025002 (1999); |hep-th/9806220| . 

U. Frisch, A. Pouquet, J. Leorat, and A. Mazure, J. Fluid Mech. 68, 769 (1975). 

M. Meneguzzi, U. Frisch, and A. Pouquet, Phys. Rev. Lett. 47, 1060 (1981). 

K. Jedamzik, V. Katalinic, and A.V. Olinto, Phys. Rev. D 57, 3264 (1998); 



|astro-ph/9606080. 



[16] K. Jedamzik, V. Katalinic, and A.V. Olinto, Phys. Rev. Lett. 85, 700 (2000); 
|astro-ph/9911100 . 



12 



[17] D. Balasara and A. Pouquet, Phys. Plasmas 6, 89 (1999). 
[18] A. Brandenburg, |astro-ph/0006186 . 



[19] H.K. Moffatt, Magnetic Field Generation in Electrically Conducting Fluids (Cambridge 
U. Press, Cambridge), 1978. 

G.B. Field, E.G. Blackman, and H. Chou, Astrophys. Journ. 513, 638 (1999); 



|astro-ph/9810345 . 



[21] K. Subramanian and J.D. Barrow, Phys. Rev. D 58, 083502 (1998); |astro-ph/9712085 . 
[22] K. Dolag, M. Bartelmann, and H. Lesch, |astro-ph/9906329| . 



13 



